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We study Casimir forces on the partition in a closed box (piston) with perfect metallic boundary 
conditions. Related closed geometries have generated interest as candidates for a repulsive force. 
By using an optical path expansion we solve exactly the case of a piston with a rectangular cross 
section, and find that the force always attracts the partition to the nearest base. For arbitrary 
cross sections, we can use an expansion for the density of states to compute the force in the limit of 
small height to width ratios. The corrections to the force between parallel plates are found to have 
interesting dependence on the shape of the cross section. Finally, for temperatures in the range 
of experimental interest we compute finite temperature corrections to the force (again assuming 
perfect boundaries). 
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I. INTRODUCTION 



A striking macroscopic manifestation of quantum elec- 
trodynamics is the attraction of neutral metals. In 1948 
Casimir predicted that such a force results from the mod- 
ification of the ground state energy of the photon field 
due to the presence of conducting boundary conditions 
The energy spectrum is modified in a fashion that 
depends on the separation between the plates, a. While 
the zero-point energy is itself infinite, its variation with a 
gives rise to a finite force. High precision measurements, 
following the pioneering work of Lamoreaux in 1997 Q, 
have renewed interest in this subject. A review of ex- 
perimental attempts to measure the force prior to 1997, 
and the many improvements since then, can be found in 
Ref. [|[. As one example, we note experiments by Mo- 
hideen et aL[J], using an atomic force microscope, which 
have confirmed Casimir's prediction from lOOnm to sev- 
eral /im, to a few percent accuracy. Forces at these scales 
are relevant to operation of micro-electromechanical sys- 
tems (MEMS) , such as the actuator constructed by Chan 
et al. [j| to control the frequency of oscillation of a nan- 
odevicc. They also appear as an undesirable background 
in precision experiments such as those that test gravity 
at the sub- millimeter scale @. 

An undesirable aspect of the Casimir attraction is that 
it can cause the collapse of a device, a phenomenon 
known as "stiction" This has motivated the search 
for circumstances where the attractive force can be re- 
duced, or even made repulsive The Casimir force, of 
course, depends sensitively on shape, as evidenced from 
comparison of known geometries from parallel plates, to 
the sphere opposite a plane [p, the cylinder opposite a 
plane [l(| , eccentric cylinde rs Till , the hyperboloid oppo- 
site a plane [ljj , a grating [l3j. a corrugated plane [l4| . 
The possibility of a repulsive Casimir force between per- 
fect metals can be traced to a computation of energy of 



a spherical shell by Boyer [15(, who found that the finite 
part of this energy is opposite in sign to that for parallel 
plates. This term can be regarded as a positive pressure 
favoring an increased radius for the sphere, if it were the 
only consequence of changing the radius. The same sign 
is obtained for a sq uare in 2-dimensions and a cube in 3 
dimensions 

[S Il7| . For a parallelepiped with a square 
base of width b and height a, the finite part of the Casimir 
energy is positive for aspect ratios of 0.408 < a/b < 3.48. 
This would again imply a repulsive force in this regime 
if there were no other energy contributions accompany- 
ing deformations at a fixed aspect ratio. Of course, it 
is impossible to change the size of a material sphere (or 
cube) without changing its surface area, and other con- 
tributions to its cohesive energy. For example, a spheri- 
cal shell cut into two equal hemispheres which are then 
separated has superficial resemblance to the Boyer calcu- 
lation. However the cut changes the geometry, and it can 
in fact be shown [l8| that the two hemispheres attract. 

The piston geometry, first considered by Cavalcanti 
frpl (in 2 dimensions) and further considered in Refs. (20l . 
[2lT | (in 3 dimensions), is closely related to the paral- 
lelepiped discussed above. 1 As depicted in Fig. [TJ we 
examine a piston of height h, with a movable partition 
at a distance a from the lower base. The simplest case 
is that of a rectangular base, but this can be generalized 
to arbitrary cross sections. This set-up is experimen- 
tally realizable, and does not require any deformations 
of the materials as the partition is moved. The force re- 
sulting from rigid displacements of this piece is perfectly 
well defined, and free from various ambiguities due to 
cut-offs and divergences that will be discussed later. In 
particular, we indeed find the finite part of the energy 
can be "repulsive" if only one of the boxes adjoining the 



1 The piston geometry was earlier mentioned in Ref. [22 



2 



partition is considered, while if both compartments are 
included, the net force on the partition is attractive (in 
the sense that it is pulled to the closest base). 

This paper expands on a previous brief publication of 
our results (2pj , and is organized as follows. Section UT1 
introduces the technical tools preliminary to the calcu- 
lations, and includes sections on cutoffs and divergences, 
the optical path approach, and on the decomposition 
of the electromagnetic (EM) field into two scalar field 
(transverse magnetic and transverse electric) with Dirich- 
let and Neumann boundary conditions (respectively). 
Details of the computation for pistons of rectangular 
cross section are presented in Section IIII| and the ori- 
gin of the cancellations leading to a net attractive force 
on the partition is discussed in some detail. Interestingly, 
it is possible to provide results that are asymptotically 
exact in the limit of small separations for cross sections of 
arbitrary shape. As discussed in Section llVl there is an 
interesting dependence on the shape in this limit, related 
to the resolution with which the cross section is viewed. 
Finally, in Scction|V]wc present new results pertaining to 
corrections to the Casimir force at finite temperatures in 
such closed geometries (for perfect metals). We conclude 
with a brief summary (Section IVIj) . and an Appendix. 

II. PRELIMINARIES 

Before embarking on the calculation of the force on 
the partition, we introduce some relevant concepts in this 
Section. Subsection III Al discusses the general structure 
of the divergences appearing in the calculation of zero- 
point energies, and indirectly justifies our focus on the 
piston geometry. The optical path approach, which is our 
computational method of choice is reviewed in Sec. Ill Bl 
Another important aspect of the piston geometry is that 
it enables the decomposition of the EM field into Dirichlet 
and Neumann scalar fields, as presented in Sec. Ill Cl 



with edges but otherwise flat, like a parallelepiped, it is 
the total length of the edges) proportional to A 3 and A 2 
respectively, and so forth. For example in the case of a 
scalar field with Dirichlet boundary conditions, we find 

E(A) = ^-1-SA^±-LA^... + E, (1) 

where denote lower order cutoff dependences, 2 and 
E is the important finite part in the limit of A — > oo. 
The EM field also enjoys a similar expansion, although 
some terms may be absent. 

Although the volume term is cancelled by an identical 
term in Eq , this is not obviously the case for the other di- 
vergent terms (surface area, perimeter, and so on). The 
energy of an isolated cavity is therefore dependent on 
the physical properties of the metal. A determination 
of the stresses in a single closed cavity requires detailed 
considerations of the metal, and its extrapolation to the 
perfectly conducting limit will be problematic [2J]. It is 
tempting to ignore these cutoff dependent terms, and to 
remove them in analogy to the renormalization of ultravi- 
olet divergences in quantum field theories. This is unjus- 
tified as there are no boundary counter-terms to cancel 
them, see Ref. (HJ. If however, we are interested in the 
force between rigid bodies, then any surface, perimeter, 
etc. terms are independent of the distance between them, 
and a well defined (finite) force exists in the perfect con- 
ductor limit. 

While the piston geometry considered in this paper is 
closely related to the parallelepiped cavities considered 
in the literature, it does not suffer from problems asso- 
ciated with changes in shape. The overall volume, sur- 
face, and perimeter contributions are all unchanged as 
the height of the partition is varied, and the force acting 
on it is finite and well defined. The same observations led 
Cavalcanti [l9| to consider a rectangular (2-dimensional) 
piston. He found that the force on the partition, though 
weakened relative to parallel lines, is attractive. 



A. Cutoff dependence 

Let us consider an empty cavity made of perfectly 
conducting material. The Casimir energy of the EM 
field is a sum over the zero point energies of all modes 
compared to the energies in the absence of the material 
E C (A) = E(A) - E (A) = ±_£ A hw m - \ £ A hw^, and 
is divergent if the upper limit A is taken to infinity. In 
a physical realization, the upper cutoff is roughly the 
plasma frequency of the metal, as it separates the modes 
that are reflected and those that are transmitted and are 
hence unaffected by the presence of the metallic bound- 
aries. Based on general results for the density of states 
in a cavity (23[, we know E(A) has an asymptotic form, 
with a leading term proportional to the volume V of the 
cavity and the fourth power of A and sub-leading terms 
proportional to its surface length L which is 

related to the average curvature of the walls (in a cavity 



B. Optical approach 

The Casimir energy can be exp ressed as a sum over 
contributions of optical paths [26|, and much intuition 
into the problem is gained by classifying the correspond- 
ing paths. For generic geometries this approach yields 
only an approximation to the exact result that ignores 
diffraction. Fortunately, it is exact for rectilinear geome- 
tries if reflections from edges and corners are properly 
included. 

Consider a free scalar field in spatial domain T> obey- 
ing some prescribed boundary conditions (Dirichlet or 



For general geometries, there are also linear and logarithmic 
terms in A, but for the class of geometries examined in this paper 
(pistons) there are no further terms in A. 
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Neumann) on the boundary B = &D. The Casimir en- 
ergy is defined as the sum over the zero point energies, 
E = ^2 \ where uj are the eigenfrequencies in D (we 
refrain from subtracting Eq for the moment). This ex- 
pression needs to be regularized, as explained in the pre- 
vious section, by some cutoff A. We choose to implement 
this by a smoothing function Sa(&) = e~ fe / A , and thus 
examine E(A) = J2 k \^{k)e~ k / K . 

The Casimir energy can be expressed in terms of the 
spectral Green's function G(x,x',fc) which satisfies the 
Helmholtz equation in T> with a point source, 



(V 2 + fc 2 )G(x', x, k) = -(5 3 (x' - x) 



(2) 



and subject to the same boundary conditions on B as 
the field. The Casimir energy of a scalar field is then 
given by the integral over space and wavenumber of the 
imaginary part of G, in the coincidence limit x' — > x [2o| 

(h = c = 1), as 



E(A) = -3* / dkk 2 e- k/A [ dxG{x,x,k). 

T Jo JV 



(3) 



The knowledge of the Helmholtz Greens function at co- 
incident points allows us to calculate the Casimir energy 
of the configuration. 

It is convenient to introduce a fictitious time t and a 
corresponding space-time propagator, G(x',x, t), defined 
as the Fourier transform of G(x',x, k). The propaga- 
tor G(t) can be expressed as the functional integral of 
a free quantum particle of mass m — 1/2 with appro- 
priate phases associated with paths that reflect off the 
boundaries. 

In the optical approach, the path integral is approxi- 
mated as a sum over classical paths of exp[iS Pr (x', x t)J , 
weighted by the Van Vleck determinant D r (x', x, t) [271 ]. 
Here S Pr (x.' , x, t) is the classical action of a path p r from 
x to x' in time t, composed of straight segments and 
undergoing r reflections at the walls. For rectilinear ge- 
ometries, like the parallelepiped that we will discuss, this 
is exact and effectively generalizes the method of images 
to the Helmholtz equation. 

For definiteness consider a scalar field satisfying either 
Dirichlet or Neumann boundary conditions, introducing 
a parameter rj, which is —1 for the Dirichlet and +1 for 
the Neumann case. The Green's function is then given 
by 



G(x',x, k) 



^ 47r/ Pr (x',x)' 



(4) 



where l Pr is the length of the path from x to x' along p r . 
There is a phase factor <fr(p r , rj) = r) ns+nc with n s and n c 
the number of surface and corner reflections, respectively. 
Note that reflections from an edge do not contribute to 
the phase. 

Since paths without reflections or with only one reflec- 
tion can have zero length, they require a frequency cut- 
off A, implemented by the smoothing function S^(k) = 



e fc / A . Then the x and k integrals can be exchanged, the 
k integral performed and the Casimir energy written as 



E V (A) 



A 4 (3 - (Z Pr (x)A) 2 ) 
(U*)A) 2 ) 3 ■ 



(5) 



The limit A — > oo can be taken in each term of the sum, 
unless a path has zero length, which can occur only for 
cases with r = or r = 1. After isolating these two 
contributions, we set 



E V {A) = E {A) + E 1 (A, ri) + E v 



(6) 



The zero reflection path has exactly zero length, and 
contributes the energy E (A) = £zVA 4 , where V is the 
volume of the space. This is a constant and therefore 
does not contribute to the Casimir force. The one reflec- 
tion paths (energy Ei(A,rf)) generate cutoff dependent 
terms, but gencrically, also cutoff-independent terms. We 
will show however that such one reflection terms do not 
contribute to the force when specialized to the piston 
geometry. 

For paths undergoing multiple reflections r > 1, the 
length l Pr is always finite, and we can safely send A — > oo 
in eq. §S§ , resulting in the simpler and cutoff independent 
contribution 



(7) 



This is a finite contribution to the energy in the limit A — > 
oo. The derivative of E^ gives the finite force between 
the rigid bodies. 



C. Electromagnetic field modes 

In the previous section we defined the optical approach 
for a scalar field. Although a similar definition can be 
made for the electromagnetic field in an arbitrary ge- 
ometry, the Helmholtz equation becomes matrix- valued, 
complicating the treatment even in a semiclassical ap- 
proximation. However, in the piston geometry, with ar- 
bitrary cross section, the EM field can be separated into 
transverse magnetic (TM) and transverse electric (TE) 
modes, that satisfy Dirichlet and Neumann boundary 
conditions. 

At the surface of an ideal conductor, the E and B fields 
satisfy the boundary conditions, E x n = and B • n = 0, 
where n is the normal vector at the surface. The normal 
modes of the piston consist of a TM set, which satisfy 



ip(y, z) cos (nirx/a) , n — 0,1,2, 



(8) 



where ip vanishes on the boundaries of the domain, 
and therefore satisfies Dirichlet conditions on the 2- 
dimcnsional boundary; and a TE set, with 



B , 



(j)(y, z) sin (mrx/a) , n = 1, 2, 3, ■ 



(9) 
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where satisfies Neumann boundary conditions. The 
other components of E and B can be computed from 
Maxwell's equations, and are easily shown to obey con- 
ducting boundary conditions. There is, however, one im- 
portant exception: the TE mode built from the trivial 
Neumann solution, <f> = constant, does not satisfy con- 
ducting boundary conditions unless the constant (and all 
components of E and B) are zero. We must ensure that 
the corresponding set of modes in eq. ([9]) are not included 
in the Casimir summation. 

Equations © and © enable us to list the spectrum of 
the electromagnetic field. Denote the spectra of the TM 
modes as the set il(Nj <g> D$) C K. Here N] indicates 
that the component on the interval satisfies Neumann 
boundary conditions, and D$ indicates that the com- 
ponent on the cross section satisfies Dirichlct boundary 
conditions. Similarly, we denote the spectra of the TE 
modes as Vt(Dj <E> N$) in a similar notation. However, as 
explained above, we must remove fl(Di), which are the 
frequencies with <p = constant. Hence, the electromag- 
netic spectra is the set 



n c = n(N T ® D s ) U f2(£> 7 ® N s ) \ il(D T ). 



(10) 



Note that fl(Dj) = {Tr/a,2ir/a, . . .} is the set of eigen- 
frequencics in 1-dimcnsion. The Dirichlct and Neumann 
spectra on the interval are identical except for the n = 
mode (see eqs. ([5} and ((9|)), but the energy of this mode is 
independent of a and does not contribute to the Casimir 
force. So we may replace iVj — ► Di in the TM spectrum 
and Dj — > Nj in the TE spectrum, with the result 

n c « n(Uj ® D s ) U tt(Ni ® Ns) \ 0(-D/), (11) 

where the notation w indicates equality up to terms in- 
dependent of a. Thus, the EM spectrum is the union 
of Dirichlet and Neumann spectra in the 3-dimcnsional 
domain, T>, except that the Dirichlet spectrum on the 
interval must be taken out. 



III. RECTANGULAR PISTON 

A. Derivation 

The piston geometry is depicted in Fig. [T] The do- 
main D consists of the whole parallelepiped, the union 
of Regions I and II. Only the partition, located a dis- 
tance a from the base and h — a from the top, is free to 
move. We study the scalar field for both Dirichlet and 
Neumann boundary conditions and the electromagnetic 
field. According to eq. (fTTj) . the EM Casimir energy arises 
from the sum of the Dirichlet and Neumann energies in 
3 dimensions minus the Dirichlet Casimir energy in one 
dimension, E = hA 2 /2tt - C(2)/(47ra) - <(2)/(4tt(/i - a)) 
(a standard result). In total, then, the EM Casimir force 
on the partition is 



Fn = F, 



D 



F, 



N 



C(2) C(2) 
Ana 2 4ir(h - a)' 




FIG. 1: (color online). The 3-dimensional piston of size h x 
b x c. A partition at height a separates it into Region I and 
Region II. A selection of representative paths are given in (a)- 
(i). Several of these paths (namely, (a,b,c,d,h,i)) have start 
and end points that actually coincide, but we have slightly 
separated them for clarity. 



where the final term vanishes if we take h — » oo. 

Let us initially focus on Region I, the parallelepiped of 
size a x b x c, below the partition. The optical energy 
receives contributions from the sum over all closed paths 
p r in domain T>i\ Each path is composed of straight seg- 
ments with equal angles of incidence and reflection when 
bouncing off the walls. There are four distinct classes 
of paths: -Ep^, from periodic orbits reflecting off faces 
(e.g. paths (c), (d), (i)); £7 aper , from aperiodic tours off 
faces (e.g. paths (a), (e), (f)); E edse , from closed paths 
involving reflections off edges (e.g. paths (b), (g)); and 
E cnx , from closed paths with reflections off corners (e.g. 
path (h)). To each path p r , we associate a vector l Pr 
pointed along the initial heading of the path, and of 
length |l p J = l Pr . 

First we consider the periodic orbits, which are paths 
that involve an even number of reflections off faces, with 
r = {0,2,4,...} (e.g. paths (c), (d), (i)). As the cen- 
tral point is varied throughout T>i, the length l Pr of 
each periodic path remains fixed, making the integra- 
tion trivial, i.e. d 3 x — > abc = V. We index the 
paths by integers n,m,l, so l Pr — (2na, 2mb, lie), with 

The n = i 



l nml = y / (2na) 2 + (2mb) 2 + (2lc) 2 . The n = m = I = 
term gives E = ^frVA 4 (see eq. (O with l Pr = 0), while 
all others are evaluated using eq. 0, giving: 



2tt 2 
3 



VA 4 



abc 
32^2 

-FA^-MA 
2ir 2 16tt 2 a 3 



Z 3 (a,b,c; 4) 

r(b,c) 



+o( 



e 2v ^ a | , as a 







(13) 



(14) 



(12) 



where Zd(a\, ■ ■ ■ , a^; s) is the Epstein zeta function de- 
fined in the Appendix (eq. ([501) 1. and T(b,c) does not 
depend on a and hence does not contribute to the force 
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on the piston. In eq. (p~4|) g = min(6, c), and A = be is the 
area of the base. The leading cutoff independent piece as 
a — > is the Casimir energy for parallel plates, coming 
from orbits that reflect off both the base and partition, 
see paths (c), (d), etc. in Fig. [I] To extract this behavior 
we have used 

Z d ( ai ,...,a d ;s) = ^- + o(-), (15) 

in the regime a± <C «2, • • • , o-d (see Appendix). 

We next consider the contribution of the aperiodic or- 
bits that involve an odd number of reflections off faces, 
with r = {1,3,5,...}. Examples in the figure include 
paths (a), (e), and (f). For each such path, when we 
vary the point of integration over Dj one of the Cartesian 
components of the length vector l Pr changes and the other 
two components are fixed. For example, only the x com- 
ponent varies for those paths that undergo an odd num- 
ber of reflections off walls parallel to the j/z-plane and an 
even number of reflections off walls parallel to both the xy 
and xz-planes. The x-component l Pr increases by 2a each 
time that the number of reflections off the yz-planes in- 
creases, so that l Pr = (2a(n— l)+2£(x), 26m, 2d), where 
£(x) = x or = a — x depending on the direction of 
the path. The summation over n and the x-integral J dx 
together combine to form an integral over x from — oo to 
+oo. So we introduce l Pr (x) = \/{2x) 2 + (26m) 2 + (2d) 2 
in terms of which the integration over the fixed compo- 
nents y and z is trivial J dydz = be, and the x-integration 
is elementary. The above example singled out the x- 
componcnt. To include all such paths, the analysis must 
be repeated for the other two components under the 
cyclic interchange of a, 6, and c. Employing eq. §5§ for 
{n, m} = {0, 0} and eq. for {n, m} ^ {0, 0} we obtain, 

#f P , r (A) = ^SA 3 --^(abZ 2 (a,b;3) + acZ 2 (a 1 c;i) 
Sir oiir \ 

+6cZ 2 (6,c;3)) (16) 

7 l n.3 C(3) P , 

= — SA^-i?— 2+T(6,c) 

Sir 647T a z 

+C) L-2ng/a\ ^ as fl ^ Q (17) 

where T does not depend on a, S — 2(ab + ac + be) = 
aP + 2A is the total surface area, and P = 26 + 2c is the 
perimeter of the base. The leading cutoff independent 
piece as a — ► comes from paths that reflect once off a 
side wall and off both the base and partition, see paths 
(c), (f), etc. in Fig. [D 

Next we calculate the contribution of even reflection 
paths which intersect an edge of Region I. Examples in- 
clude (b) and (g) in Fig. [TJ In this case it is only the 
component of l Pr parallel to the edge that remains fixed, 
while the other 2 components vary as the point of origin 
varies over T>i. For example, suppose the reflecting edge 
is parallel to the z-axis. Then the z-integration is trivial, 
J Q C dz — c, and the path vector is a function of x and y 
given by l Pr = (2a(n-l)+2£(x), 2b(m-l) + 2ip(y), 2d), 



where £(x) = x or £(x) = a — x, ip(y) = y or ip(y) = b — y, 
depending on the quadrant that l Pr lies in: up or down 
in x, right or left in y, respectively. In this case we can 
replace both the summations over n and m and the dou- 
ble integral over x and y by an integral over the whole 
xy-planc. This integral is most easily performed in polar 
co-ordinates, using a path length that may be written 
as l Pr (r) = \J (2an) 2 + (2r) 2 . The contribution to the 
Casimir energy is found to be 

£UA) = J_^-f)(I + UI), as, 

g 327T lo7r \a b c J 

where L — 4 (a + 6 + c) = 4a + 2P is the total perimeter 
length. The cutoff independent piece ~ 1/a comes from 
paths that reflect once off a side edge and off both the 
base and partition, see path (g), etc. in Fig.[TJ 

Finally, we consider the paths which reflect off a cor- 
ner (E CIlI ). In this case, as the integration variable moves 
throughout its domain, all components of the distance 
vector l Pr vary. Hence, we can incorporate all such paths 
by extending our integral over all space in x, y and z. 
This leaves no dependence on the geometry of the paral- 
lelepiped (i.e., it is independent of a, 6, and c), and only 
contributes a constant that is of no interest, which we 
ignore. 

Adding together all these contributions, we obtain the 
Casimir energy of a scalar field in Region I as 

E^A) = Ay A 4 + AsA 3 + -LiA 2 + E*, (19) 

where E^ gives the cutoff independent piece, from 
eqs. ([13]), (HU), and fig]). We note that the cutoff de- 
pendent terms agree with the leading terms obtained by 
integrating Balian and Bloch's asymptotic expansion of 
the density of states [23|. 

We obtain the Casimir energy for the entire piston by 
adding to eq. fTi?]) the analogous expression for Region 
II obtained by the replacements: a — > h — a, V — > hA, 
S — > hP + 4A, and L — > 4h + 4P. It is easy to see that 
after including Region II the sum of all cutoff dependent 
terms is independent of partition height a. Therefore 
the force on the partition is well defined and finite in 
the limit A — > oo. Also, of course, the contribution to 
the Casimir energy from the region outside the piston is 
independent of a and can be ignored entirely. The force 
on the partition is given by the partial derivative with 
respect to a of the cutoff independent terms as 

F v = —^ b, c) + E v (h - a, 6, c)) , (20) 

where we have defined E^ = E v (a, 6, c). 

We focus on the h — > oo limit in which the expression 
for the contribution from Region II simplifies. Consider 
the periodic, aperiodic, and edge paths whose cutoff inde- 
pendent contribution to the energy is given in eqs. Q13p . 
(HHJ), and (fT5|) . Replacing a — > h— a, taking h — > oo, and 
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using cq. (|52|) of the appendix in these equations gives 
h — a 



E 



n 



E 



11 



h - a ( 1 1\ , , 

+ ? c(3) ' 



0, 



(21) 



where we have not reported terms independent of a, since 
they do not affect the force. Also, we re-express the Re- 
gion I energy E v (a,b,c) in a fashion that is useful for 
a <C 6, c, using cq. (|51[) of the appendix. The net force 
on the partition due to quantum fluctuations of the scalar 
field is then 

_3C(4) A C(3) P C(2) J,(6/c) 



16tt 2 a 4 



32tt a 3 167ra 2 32tt 2 A 



77^-^ ^ n2 (Ko (2irmnb/a) b + (b <-> c)) 



m,n— 1 



coth(/ mn (6/a, c/a)) 



32 a " ^ /mn(&/a, c/a) sinh 2 (/ mn (6/a, c/a)) 



,(22) 



where the primed summation is over {m, n} G Z 2 \ {0, 0} 
and Kq is the zeroth order modified Bessel function 
of the second kind. Here we have defined J v (x) = 
Z 2 ( x 1 / 2 1 x- 1 ' 2 -A) + nrjjx + .t' 1 )^) and f mn (x,y) = 
TT^J(mx) 2 + {ny) 2 . The first four terms of eq. ([22]) dom- 
inate for a <C b,c, while the following terms are ex- 
ponentially small in this regime. The first term arises 
from periodic orbits reflecting off walls (see cq. fll4[l ). 
the second term from aperiodic tours bouncing off walls 
(see eq. (|17|l). the third term from reflections off edges 
(see eq. the fourth term from Region II paths (see 

eq. (f2"Tj) ). Note that the infinite series, involving expo- 
nentially small terms, is convergent for any a, b, c. 

The electromagnetic case is closely related to the scalar 
Dirichlet and Neumann cases, which we discussed in de- 
tail in Section III CI According to cq. (TT2"j). the EM 
Casimir energy in Region I is related to the Dirichlet 
energy Ed and the Neumann energy En by 

E I C (A) = E I D (A) + E I N (A)- ]T E 1 (d,A), (23) 

d—a,b,c 

where E\(d,A) = dA 2 /2n - (^(2) /(And) is the energy of a 
scalar field in 1-dimensions obeying Dirichlet boundary 
conditions in a region of length d. The contribution from 
Region II follows from replacing a — > h — a. Combin- 
ing previous results, the electromagnetic Casimir force is 
found to be, 
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3C(4)A { C(2) Jc(b/c) 
8tt 2 

A. v 



32tt 2 ^ 



7T" 

Tg 



coth(/ mn (6/a,c/a)) 



a m~ti fmn(b/a, c/a) sinh (f mn (b/a,c/d)) 



,(24) 



where J c (x) = J_ x (x) + J+ 1(1) = 2Z 2 (x 1 / 2 , a;" 1 / 2 ; 4). 




FIG. 2: (color online). The force F on a square piston 
(b = c) due to quantum fluctuations of a field subject to 
Dirichlet, Neumann, or conducting boundary conditions, as 
a function of a/6, rescaled as F' = 16tt 2 AF/(3C(4)) (F' = 
8?r 2 J 4F/(3C(4))) for scalar (EM) fields. The solid lines are for 
the piston, solid middle = Fc, solid upper = Fn, and solid 
lower = Fn, while their dashed counterparts are for the box. 



B. Discussion 

Here we address the implications of eqs. (f2"2"j) and 
for the force on the partition in more detail. To begin, 
we discuss the important issue of attraction versus re- 
pulsion. We are interested in comparing the force on the 
partition (Fr, where T — D, N or C for Dirichlet, Neu- 
mann, or EM boundary conditions respectively) to the 
force reported in the literature for a single cavity, which 
we denote -Fi^box [13 ■ The latter is obtained by the fol- 
lowing prescription: calculate the energy in a single rec- 
tilinear cavity, drop the cutoff dependent ( "divergent" ) 
terms, ignore contributions from the region exterior to 
the cavity, and differentiate with respect to a to obtain 
a force. We emphasize that there is no justification for 
dropping the cutoff dependent terms, so although we re- 
fer to this result, for convenience, as -Fbox, it does not 
apply to the physical case of a rectilinear box. 

For the piston geometry, we note that the sole contri- 
bution from Region II is the a-independent term denoted 
by J. In fact this is the only term that distinguishes F 
from F hox , i.e., 



F r = F r ,box - J r (6/c)/(32^4). 



(25) 



Naively, the difference by a constant may not seem impor- 
tant. Indeed it is not too important for small values of the 
ratio a/ (b, c). However it is very important for a > (b, c). 
In Fig. [5] we plot both these forces for a square cross sec- 
tion (b = c) as a function of a/6. (The plots include scalar 
as well as EM cases.) Note that in all cases F — * 0, while 



F h 



J(l)/(32n 2 A) (a constant) as a/6 



For 



this geometry J D (1) w -1.5259, J N (l) w 13.579, and 
Jc(l) ~ 12.053, so J is negative for Dirichlet and posi- 



tive for both Neumann and EM. We see that F is always F/F\\ 



attractive, while ft, ox can change sign. It is always at- 
tractive for Dirichlet, but becomes repulsive for Neumann 
when a/b > 1.745 and for EM when a/b > 0.785. This is 
the consequence of ignoring Region II and the cutoff de- 
pendence. Indeed, it is easy to show that the piston force 
is attractive for any choice of a, 6, c, h. A final comment 
is that for h finite and a = h/2, the partition sits at an 
unstable equilibrium position. This comment was made 
in Ref. [28| , although the above detailed results were not 
derived there. 

With the explicit form for F, we can more closely 
compare the piston with Casimir's original parallel plate 
geometry. For better comparison in Figs. [3] and 31 wc 
have plotted the forces for the scalar and EM fields, 
after dividing by the parallel plates results, -f^z^jv)!! = 
-3C(4)A/(167r 2 a 4 ) or F c \\ = -3C(4)A/(87r 2 a 4 ). ' First, 
note that for the EM case not only does Fc — > as 
a/b — > oo but it does so rather quickly. Since Fqu van- 
ishes as 1/a 4 , it is clear from Fig.[4]that Fc vanishes even 
more rapidly. In fact it vanishes exponentially fast, as 
e -2ira/b f or a ^ 5 We can understand this as follows: In 
this limit the most important paths are those that reflect 
off the top and bottom plates once, and therefore travel 
a distance 2a. The transverse wavenumber k = 7r/& due 
to the finite cross section, acts as an effective mass for 
the system, and damps the contribution of these paths 
exponentially. In fact for any rectangular cross section 
wc find 



Fr 



1 



_ e -2wa/b _|_ 



1 



3 — 2-Ka/c 



as a 



(26) 

Experimentally, values of a/b ~ 1 are not yet realiz- 
able. Instead, typical experimental studies of Casimir 
forces have transverse dimensions that are roughly 100 
times the separation between the "plates". This means 
that the leading order corrections to Fq\\ are more likely 
to be detected experimentally. In Figs. [3] and [4] we 
show the result of including successive corrections to 
Fu for scalar and EM cases, respectively; we plot the 
curve which includes {1/a 4 , 1/a 3 , 1/a 2 } terms and an- 
other curve that includes {1/a , 1/a 3 , 1/a 2 , 1} terms. In 
the EM case we note that the 1 /a 3 term that appears in 
the expansion for Dirichlet and Neumann boundary con- 
ditions is canceled. In the next Section, we will demon- 
strate that this is a general phenomenon for any cross 
section (see ahead to eq. (J2ZJ))- Hence the first correc- 
tion to the EM result scales as 1/a 2 , which is 0(a 2 /A) 
compared to -Fbll- We see that this correction is quite 
accurate up to a/b ~ 0.3. We suspect this regime of 
accuracy to be roughly valid for any cross section. 
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FIG. 3: (color online). The force F on a square partition (b = 
c) due to quantum fluctuations of a scalar field as a function 
of a/b, normalized to the parallel plates force F\\. Left figure 
is Dirichlet; solid middle = Fn (piston), dashed = Fo,ho* 
(box), solid upper = {1/a 4 , 1/a 3 , 1/a 2 } terms, solid lower = 
{1/a 4 , 1/a 3 , 1/a 2 , 1} terms. Right figure is Neumann; solid 
middle = JW (piston), dashed = -Fjv, b ox (box), solid lower = 



{1/a 4 
terms. 



1/a 3 , 1/a 2 } terms, solid upper = {1/a 4 , 1/a 3 , 1/a 2 , 1} 



an asymptotic expansion for small separation a. The 
generalized piston maintains symmetry along the vertical 
axis, and its geometry is the product of I ® <S of the 
interval / = [0, a] and some 2-dimensional cross section 
5 c M 2 . Let us denote by £ = k 2 the eigenvalues of 
the Laplacian on the piston base S and the interval I, 
separately, with appropriate boundary conditions 



IV. GENERAL CROSS SECTIONS 

The piston for arbitrary cross section cannot be solved 
exactly, but we can obtain much useful information from 



I- 



(27) 



The corresponding densities of states are denoted by by 
pi and ps, respectively. Then the density of states (per 
unit "energy", £) of the problem in the 3-dimcnsional 
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FIG. 4: (color online). The force F on a square partition 
(b — c) due to quantum fluctuations of the EM field as a 
function of a/6, normalized to the parallel plates force Fj|. 
Solid middle = Fc (piston), dashed = Fc.hax (box), solid lower 
= {1/a 4 , 1/a 2 } terms, solid upper = {1/a 4 , 1/a 2 , 1} terms. 



region I®S is p{£ ) and can be written as the convolution 



d£' Ps {£-£')pi{£'). 



(28) 



The 2-dimcnsional density ps is not known in general, 
since the wave equation can not be solved in full general- 
ity in an arbitrary domain S. However, for small height 
to width ratios, the smallness of a translates to high en- 
ergies £, and we will see that the asymptotic behavior 
of ps is sufficient for extracting an asymptotic expansion 
for the force is powers of 1/a. 

The number of eigenstates with energy less than £ in 
S has the asymptotic expansion at large £ [29[, 



N s (£) 



^ + V^V£ + x + r N (£)^e(£). 



(29) 



Here, x 1S related to the shape of the domain S through 



^24 (a; ir) + ^ I2ir J 7] K 



( 7j -)d7j, (30) 



where is the interior angle of each sharp corner and 
ft(7j) is the curvature of each smooth section. It is easy 
to check that \ = 1/4 for a rectangle and \ — 1/6 for 
any smooth shape (for example a circle). Note that we 
have included the step function 0(£) in the expression 
for Ns(£), ensuring that only £ > contributes. Here 
rjsr(£) is a function which designates lower order terms 
(remainder) in an £ — ► 00 asymptotic expansion. For 
any polygonal shape r/v is exponentially small, r^{£) = 
0(e~ c£ ) (c > is a constant) [30] ■ However, we are aware 
of only a much weaker estimate for smooth shapes, as 



r N {£) = <D(1/V£) H 
density of states 3 



The derivative of Ns(£) is the 



.4 
4tt 



where we have used &'(£) = S(£), and £ S(£ ) = 
V£5(£) = for all£. 

The other function in the convolution, the 1- 
dimensional density of states, is known exactly: it is sim- 
ply a sum over delta functions, which we rewrite in terms 
of its Poisson summation 



PI (£) = J2S[£ 

n=l 
a <d(£) 



2n V£ 



2 / dx cos(2irmx) S { £ 



(32) 
.(33) 



The hrst term in eq. (|33[) is the smooth contribution to 
the density of states, and the second is the oscillatory 
component. The leading contributions (as a — > 0) to the 
3-dimensional density of states come from convolving the 
smooth part of p$ with pi, giving^ 



a (-^AV£ + ^j-P + + n(£) 

\4tt 2 16tt 2tt y/s 



OO 

E 

m— 1 

ax 
ttV£ 



( — 4— sm(2maV£) -+ 
008(277^0-^) + t%{£) 



aP 

rj— — Jo(2mav / f ) 



8tt 



(34) 



The hrst line agrees precisely with the Balian and Bloch 
theory of the density of states [23| , and gives the cutoff 
dependent terms in the Casimir energy 



E(A) 



1 f°° 

=- / d£p{£)V£e-^^ A . 

2 Jo 



(35) 



The cutoff dependent contributions have no effect on 
the Casimir force when Region II is included, since they 
are linear in a, as explained earlier. The second line in 
cq. (|34p gives the leading three terms in an asymptotic 
expansion of the force 



' 167r 2 a 4 



V 



C(3) D C(2)X 



327ra 3 



;P - 



4:ira 2 



+ r v (a). (36) 



Also, even for these general cross sections, the EM energy 
can be related using eq. (fT2"|) to Dirichlct and Neumann 
energies, as 



rc—m* 

87r 2 a 4 



C(2)(1-2 X ) 
Ana 2 



rc{a). 



(37) 



3 In Eqs. d 3 1 1 ) and I I34I I we have denoted the various remainders 
by r p (£), ri(£), ri(£. ). We discuss the size of the remainders in 
the expansion of the forces F v & Fq following eq. H37I I. 
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In eqs. ([36]) and (|37|) we have written the remainder terms 
as r v (a) and re (a) (= r_i(a)+r +1 (a)), respectively. Fol- 
lowing our earlier estimates for rjy(£) that appears in 
Ns(£), and noting that there is always an 0(1) term that 
comes from Region II, we have r Vt c(a) = 0(1) for polyg- 
onal shapes and r^^cio) = 0(l/a) for smooth shapes, as 



The generalization to arbitrary cross sections in 
cq. (|37p has interesting features. The correction to the 
parallel plates result depends on geometry through the 
parameter x, which depends sensitively on whether the 
cross section is smooth or has sharp corners. For ex- 
ample, x = | for all smooth shapes and x = f° r 
an n-sided polygon of equal interior angles. Given un- 
avoidable imperfections in any experimental realization, 
one may wonder what precisely constitutes "smooth" or 
"sharp." Note that for any deformation with local ra- 
dius of curvature R (R = for perfectly sharp corners), 
we have the dimcnsionless quantity R/a, where a is the 
base-partition height. Given that our expansion is valid 
for small a, we conclude that R 3> a is a smooth deforma- 
tion, while R <^ a can be regarded as a sharp corner. As 
a simple example, consider a shape that is roughly square 
(4-sidcd polygon) if viewed from large distances, but is in 
fact rounded with radius R at the "corners" if examined 
closely Let us also imagine that the overall width (b) is 
much larger than R. Then, since the corresponding term 
in the Casimir force goes as 1 — 2\ (see eq. (|3"T)) ), we ex- 
pect the correction to be £(2)/(67ra 2 ) for a/R <C 1 and 
decrease to £(2)/(87ra 2 ) for a/R 3> 1. A more interest- 
ing example would be a self-similar (fractal or self-affinc) 
perimeter, in which the number of sharp corners deceases 
as a power of the resolution a. For such a case, we expect 
a correction scaling as a non-trivial power of 1/a, remi- 
niscent of results in Ref. [3l[ . It would be interesting to 
see if such corrections are experimentally accessible. 



Another noteworthy feature of eq. ([37)) is that the 
leading correction to the EM force (compared to par- 
allel plates) is smaller by order of a 2 /A. By contrast the 
corrections are only of order aj\J~A for scalar fields with 
either Dirichlet or Neumann boundary conditions. How- 
ever, the latter corrections are exactly equal and opposite 
in sign, and cancel for the EM force. It is interesting to 
inquire if this precise cancellation applies only to perfect 
metallic boundary conditions, or remains when the ef- 
fects of finite conductivity are taken into account. More 
work is necessary to understand the finite conducting pis- 
ton. Yet another case is for side walls made of dielectrics, 
where a simple modification of the optical path method, 
which replaces the sign factor r\ with the reflection coef- 
ficients for TM and TE modes, suggests that the cancel- 
lation does not occur. A piston that is made entirely of 
a uniform dielectric is examined in Ref. 321 



V. THERMAL CORRECTIONS 

The question of the leading corrections to the Casimir 
force at finite temperatures T has generated recent inter- 
est, both from the practical need to evaluate the accu- 
racy of experiments, and due to fundamental issues. In 
particular, there is controversy pertaining to the appro- 
priate model for the metallic walls, which we shall ignore 
in this chapter. Instead, we shall compute corrections to 
the Casimir force due to finite temperature excitations 
of the modes in the piston, while continuing to treat its 
walls as perfect metals [33j . 

A. Rectangular piston 

We first answer this question for the piston with rect- 
angular cross section. In units with H = c = ks = 1 ; the 
inverse temperature (3 = l/T introduces a new length 
scale whose size relative to the dimensions a, 6, and 
c of the piston (we imagine, as earlier, that h — > oo) 
sets the importance of thermal corrections. (More pre- 
cisely, 7r/3 is the appropriate length scale.) In typical 
experiments a ~ Ifim, b,c ~ 100/wn, and at room 
temperature 7r/3 ~ 20fj.m. Thus the regime of most 
experimental interest is where the length scales satisfy 
a <C tt(3 < b, c. In light of this we focus on thermal 
lengths much larger than the base-partition height, i.e. 
a <C it/3. To fully investigate the low temperature regime, 
we assume a <C tv/3, b, c <C h, but will allow 7r/3 to be less 
than or greater than b or c. 

Each mode of the field can be regarded as an inde- 
pendent harmonic oscillator, and by summing the corre- 
sponding contributions, we find the free energy 

T tot = — £ In ^ _;_ gMm j =E + 6F. (38) 

We have separated out the the zero-temperature Casimir 
energy E, from the finite temperature "correction" ST = 
SE — T5S, and focus on the latter for calculating finite 
temperature effects. 

First, a note of caution is in order regarding the scalar 
field with Neumann boundary conditions. In any cav- 
ity, there is a trivial solution to the Neumann problem, 
namely a constant field with u> = 0. This means that 
whenever (3 is finite (T > 0) then ST = — oo, which 
signals condensation of the scalar field into the ground 
state. We note that this phenomenon occurs for closed 
geometries where the spectrum is discrete and not in gen- 
eral for open geometries in which the region near u> = 
is integrable due to phase space suppression. We will 
proceed by calculating the free energy of a scalar field 
with both Dirichlet and Neumann boundary conditions, 
ignoring the mode with oj = for the latter. We then 
use eq. (|T2|) to obtain the EM force. This procedure is 
valid since the offending Neumann mode is specifically 
excluded from the EM spectrum. 
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For a Dirichlct scalar field in Region I. since all modes 
satisfy u) m > ir/a, their Boltzmann weights are small in 
the limit of a <C 7T/3, and 



5T 1 = O 



--rr/3/a 



(39) 



is exponentially small. Similarly the a-dependent terms 
of the electromagnetic free energy in region I are expo- 
nentially small. This is true for any cross section and 
reflects the fact that thermal wavelengths ~ tt(3 are ex- 
cluded from Region I (34|. However, a power law contri- 
bution to the free energy and force will come from Region 
II. We use the optical expansion, which remains exact for 
the free energy in rectilinear geo metries, to compute this 
contribution for scalar fields [34| , as 



q 



[U*) 2 



1 



(<7/3) 5 



(40) 

Note that here the sum ranges over 5 e Z \ {0} — the 
q = term is just the Casimir energy (see eq. (J7])). 

It is natural to break the energy up into the famil- 
iar four classes of paths: periodic orbits, aperiodic tours 
off faces, reflections off edges, and reflections off corners. 
However, summing each set separately gives a logarith- 
mic divergence (that cancels among the different classes 
for Dirichlet boundary conditions ). Fortunately, this 
problem can be ignored in the h — ► 00 limit, as can be 
seen, for example, by considering the contribution from 
the sum over periodic orbits (paths (c), (d), etc in Fig. 1). 
Noting that ft, — a is the height of the piston in Region 
II, we have 



00 00 



ST 
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W2 E E 

q—1 n,m,l— — oo 

(h — a)bc 



16tt 



[(n(h - a)) 2 + (mb) 2 + (lc) 2 + {qf3/2f 



.(41) 



This expression is logarithmically divergent, but if we 
take h — > 00, only the n = term contributes and the 
remaining summation over {q, m, 1} is finite. Strictly 
speaking, the interchange of the limit h — > 00 with 
the summations, which eliminates the logarithmic diver- 
gence, is formally problematic. However a more rigorous 
analysis justifies this step for the Dirichlet case through 
the cancellation among the different classes, but always 
gives —00 for the Neumann case as anticipated. Per- 
forming this interchange gives the following result for the 
contribution of periodic orbits 

xr ii gmVp-Aa) {±-a)A 



(h - a) A 
32tt 2 ' 



,(&,c;4) + 0(V 47r9/ ' 3 ) (42) 



with g = min(&, c) and V p as the total piston vol- 
ume. Here we have expanded for small (3 relative to 



g = min(6, c). We note that although the third term 
is independent of (3, this really is part of 6 J 7 . The reader 
that is interested in the opposite limit of (3 — > 00, i.e., 
the low temperature limit, should look ahead to Section 

ESI 

Proceeding in a similar fashion for all contributions 
to the free energy of a scalar field we find (ignoring the 
exponentially small contribution of Region I) 



5T n 



((4)(V p ~Aa) {(3)(S p ~Pa) 

n 2 (3 4 11 8tt/3 3 

C(2)(/i-o) M n (b/c)(h-a) 



4tt/3 2 32tt(3VA 
J 71 {b/c){h~ a) TT 2 (V P -Aa) 



32ir 2 A 



8/3 4 



^ , 1 + 2f mn (b, c) - e~f"^c) (j, _ a) 



- n /„3„(&,c)sinh 2 (/ m „(6,c)) ■' P 2 

OO 

^2 — (K^irmnb) + Ki(2irmnc)) (43) 



where we have defined M n (x) 



1/2 



-1/2 



Z 2 (x 3 / 2 ,x-V 2 ;3) 



)C(2), b = 2b/f3,c = 2c/ (3, and S p is the 



total surface area of the piston. It is important to note 
that while <5J r _ 1 = SJ-fj, 5T + \ = 5Tn is not strictly 
correct as we have ignored the ui = Neumann mode. 
Although 5Tn = —00, as stated earlier, this expression 
correctly gives the a-dependencc in STn- 

The EM case can be handled in a similar fashion. Re- 
peating our earlier decomposition, we note that SJ- EM = 
8T-1 + 8T + i + ((2)(h — a)/f3 2 , since the spectral decom- 
position in eq. (fTTj) correctly leaves out the w = mode 
of the Neumann spectrum. We thus find (again ignoring 
the exponentially small contribution of Region I) 



5T E 



2C(4)(V„- Aa) C(2)(h-a) 



ir 2 /3 4 
M c (b/c)(h-a) 



2tt(3 2 
J c (b/c)(h-a) 



32Tr(3y/A Z2ir 2 A 
tt 2 (V p -Aa) s ^,l + 2f mn (b, c) - e - 2 />"-(M 



4/3 4 



E' 



/mn(^c)sinh (f mn (b,c)) 



(44) 



where M c (x) = M_i(x) + M +1 (x) = 2Z 2 (x 3 / 2 , x~ 3 / 2 ; 3). 

In Eqs. (H2]) and (|44p we have written the expansion 
as a series in increasing powers of (3. The result, though, 
is correct (up to exponentially small terms in a/ir(3) for 
any ratio of (3 to b or c, and for h much larger than any 
of the other scales. The infinite summations that appear 
are convergent for all finite values of {/?, 6, c}. The lead- 
ing term in eq. (|44|) is the Stefan-Boltzmann energy, and 
the following terms are corrections due to geometry The 
term independent of (3 is equal to but opposite in sign to 
that appearing in the Casimir energy. Note that the ap- 
pearance of a term independent of (3 is an artifact of per- 
forming a small (3 expansion. All terms depend linearly 
on a and provide a constant force on the partition. Note 
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that the first five terms in ST^ and the first four terms 
in ST EM have power law dependences on /?, while the re- 
maining terms (summations) are exponentially small for 
7T/3 < (b, c). 



B. General cross section 

If we consider general I <E) S geometries, as in Sec- 
tion IIV1 we may use the smooth 3-dimcnsional Balian 
and Bloch density of states in Region II to obtain the 
leading terms in the free energy. Specifically, we use the 
first line of eq. (|34|) with the replacement a — > h — a for 
p{£), and calculate the free energy from 



f t /f> 



- / d£p(8)]n(l-exp(-pV£)). 
o Jo 



(45) 



Since we only know the first three terms in the expan- 
sion for the density of states, we will obtain contributions 
proportional to the volume, surface, and perimeter of the 
piston, but nothing at 0(1/ (3). It is fairly straightforward 
to get 



5T L 



C(4)(V p -Aa) C(S}{Sp-Pa) 



7T 2 /3 4 



a) 



TTf3 2 

2((4:)(V P -Aa) 



7T 2 (3 4 
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+ 



87T/3 3 

1" 

C(2)(l-2 X )(ft-o) 



(46) 



7T/3 2 



-a 



(47) 



We again see the effect of the modes excluded from Re- 
gion I due to a <C n(3, in the factors V p — Aa, S p — Pa, 
and h — a. These leading terms provide thermal contri- 
butions to the quantum force on the partition, given in 
eqs. ([go]) and ([57|1. 

Let us comment on the relationship between the 
Casimir and thermal contributions to the force. We be- 
gin by focusing on the regime that is perhaps of most 
experimental interest: a -C 7T/3 <C \[A. If we include 
both Casimir and thermal contributions to the force, as 
given in eqs. (J37J) & (|47)) . 



3C(4) / 1 



1 



8tt 2 \a 4 (f3/2) 4 
C(2)(1-2 X ) ( 1 

47T 



1 



W2)V - (48) 

Note that the leading contributions are related to terms 
in the Casimir energy by the interchange a <-> (3/2, but 
this connection ceases for higher order corrections. We 
have only calculated further terms for the parallelepiped 
and we can compare them in this limit. In particular, 
eq. ([4*4")) includes a contribution of order 1/(3 which has 
no counterpart (i.e. a term of order 1/a) in the Casimir 
force. A term of order 1/a can only come from the deriva- 
tive of ~ In a, which is absent from the EM Casimir en- 
ergy. 
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FIG. 5: (color online). The force Ft from thermal fluc- 
tuations on a square partition (b — c), normalized to 
the Stefan-Boltzmann expression Fsb = — C(4)A/(7r 2 /3 4 ) 
(-2C,(4)A/(-k 2 (3 4 )) for Dirichlet (EM) fields, as a function of 
(3 lb. This is valid in the regime: a {n(3, b, c}. Starting from 
a normalized value of 1, the full result for Dirichlet (electro- 
magnetic) is the lower (upper) curve. Also, starting from a 
normalized value of 0, the exponentially small asymptote (as 
(3/b — » cxj) for Dirichlet (electromagnetic) is the lower (upper) 
curve. 



C. Low temperature limit 



Equation (|48j) provides the leading terms in the 
Casimir force in the limit a <C 7r/3 <C b, c (or more gener- 
ally a -C 7T/3 <g; \f~A for non-rectangular cross sections). 
We may more accurately refer to this as a "medium tem- 
perature" regime, as opposite to a lower temperature 
regime with 7r/3 ^> \J~A. In fact, for the rectangular pis- 
ton we obtained in eqs. (|43|) and (|44|) results that are 
valid for a <C \k(3, b, c} for any ratio of (3 to b or c, and 
will now examine their lower temperature limit. A naive 
application of the proximity-force approximation gives 
always a thermal correction to the force that vanishes 
as ~ 1/J3 4 = T 4 in the T -> limit [l6[. However, in 
Ref. [34[ it is argued that for open geometries this limit 
is quite subtle and is sensitive to the detailed shape of 
each surface. In fact it is reasonable to argue that for 
the cases relevant to experiments there may be weaker 
power laws, i.e., l/(3 a with a < 4. But in our closed 
geometry another scenario is natural: If T —> 0, so that 
(3 3> a, \f~A, modes are excluded from both regions due 
to a gap in the spectrum, resulting in an exponentially 
small free energy, which (for the rectangular piston) is 

V2(3 3 / 2 \Vb sfc J 

as (3 — ► oo. (49) 

A plot of the force Fxr = —dSJ-r/da (where r = D or 
C as for T = 0), derived from eqs. (|4"3"|) and ([i4]) is given in 
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Fig. [5] The force is normalized to the Stefan-Boltzmann 
term, F SB = -((4)A/ (tt 2 (3 A ) (-2C(4)A/(tt 2 /? 4 )) for 
Dirichlet (EM) fields. Having taken a <C 7T/3 in our anal- 
ysis, the a dependence is ignorable, and we plot the force 
as a function of (3/b (b = c). The high (3/b asymptotic 
curves (eq. (|49| is the EM case) are also included. Note 
that from eq. ((35]) we can read off the small P/y/A cor- 
rections to Fsb for arbitrary cross sections. 
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VI. CONCLUSIONS 



APPENDIX 



In this work we have obtained an exact, analytic re- 
sult for the Casimir force for a piston geometry Exact, 
analytic results are rare in this field but nonetheless par- 
ticularly useful for comparison with the approximations 
needed to describe real systems and more complicated 
geometries. 

We have obtained analytic expressions for the force 
acting on the partition in a piston with perfect metal- 
lic boundaries. The results are exact for the rectangular 
piston, and in the form of an asymptotic series in 1/a 
for arbitrary cross section. We find that the partition is 
always attracted to the (closer) base: consistent with a 
more general result obtained in Ref. Since the pis- 
ton geometry is closely related to single cavity for which a 
repulsive force has been conjectured, we are able to shed 
some light on this question. In particular, we empha- 
size that to avoid unphysical deformations (and closely 
related issues on cutoffs and divergences) it is essential 
to examine contributions to the force from both sides of 
the partition. The cutoff independent contribution from 
a single cavity (that we call F box ) approaches a constant 
for large a. However, in the piston geometry compensat- 
ing contributions from the second cavity cancel both the 
cutoff dependent terms and part of the cutoff indepen- 
dent term, to cause a net attraction. 

For general cross sections we find interesting depen- 
dence on geometrical features of the shape, such as its 
sharp corners and curved segments. We have obtained 
the first three terms for scalar fields and the first two 
terms for EM fields (one less due to cancellation) in an 
expansion in powers of a. By comparison to our cal- 
culated exact result for a rectangular cross section we 
estimate that this expansion is valid for a/6 « 0.3. This 
covers the conventional experimentally accessible regime, 
and is therefore a useful result for a large class of geome- 
tries. We have also obtained thermal corrections which 
cover the experimentally accessible regime. 



The general Epstein Zeta function is defined as 
Z d (ax, . . . ,a d ;s) = ^ ' ((nidi) 2 H h (n d a d ) 2 ) 



m,...,nd 



(50) 

where the summation is over {rii, ...,n d } G Z d \ 
{0, ...,0}. Note that the Riemann Zeta function is a 
special case of this, namely £(s) = Z\{\\s)j e X. 

In eq. (fTo) we pointed out that such functions could 
be approximated by a term involving the Riemann zeta 
function and a power of a±, as a± — > 0. An exact repre- 
sentation, as discussed in Ref. [l7[, is 



Z d (ai, ...,a d ;s) 

_ 2C00 , r(*=*)0F 



r(f) 0l 



Z d -i(a 2 , ...,a d ;s-l) 



At: 8 ' 2 



)/2 



x K (s-l)/2 2tto 



V( fl 2«2) 



+ {a d n d )' 2 



a 1 



\J{a 2 n 2 y H h (a d n d y 



(l-s)/2 



01 



(51) 



where K v is the modified Besscl function of the second 
kind. This is useful in Region I where a,\ is small (with 
ai — * a). 

For Region II it is important to examine the limit in 
which one of the lengths is infinite, say a\ — > 00 (with 
a,\ — > h — a). In this limit the order of the zeta function 
is reduced: 

Z d (ai, . . . , a d ; s) — ► Z d _i(a 2 , . . . , a d ; s). (52) 
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